Cohort profile: BioMD-Y (biopsychosocial factors of major depression in youth) - a biobank study on the molecular genetics and environmental factors of depression in children and adolescents in Munich.

PURPOSE: BioMD-Y is a comprehensive biobank study of children and adolescents with major depression (MD) and their healthy peers in Germany, collecting a host of both biological and psychosocial information from the participants and their parents with the aim of exploring genetic and environmental risk and protective factors for MD in children and adolescents. PARTICIPANTS: Children and adolescents aged 8-18 years are recruited to either the clinical case group (MD, diagnosis of MD disorder) or the typically developing control group (absence of any psychiatric condition). FINDINGS TO DATE: To date, four publications on both genetic and environmental risk and resilience factors (including FKBP5, glucocorticoid receptor activation, polygenic risk scores, psychosocial and sociodemographic risk and resilience factors) have been published based on the BioMD-Y sample. FUTURE PLANS: Data collection is currently scheduled to continue into 2026. Research questions will be further addressed using available measures.

Human cortical neurogenesis is altered via glucocorticoid-mediated regulation of ZBTB16 expression.

Glucocorticoids are important for proper organ maturation, and their levels are tightly regulated during development. Here, we use human cerebral organoids and mice to study the cell-type-specific effects of glucocorticoids on neurogenesis. We show that glucocorticoids increase a specific type of basal progenitors (co-expressing PAX6 and EOMES) that has been shown to contribute to cortical expansion in gyrified species. This effect is mediated via the transcription factor ZBTB16 and leads to increased production of neurons. A phenome-wide Mendelian randomization analysis of an enhancer variant that moderates glucocorticoid-induced ZBTB16 levels reveals causal relationships with higher educational attainment and altered brain structure. The relationship with postnatal cognition is also supported by data from a prospective pregnancy cohort study. This work provides a cellular and molecular pathway for the effects of glucocorticoids on human neurogenesis that relates to lasting postnatal phenotypes.

Methylation Patterns of the FKBP5 Gene in Association with Childhood Maltreatment and Depressive Disorders.

Childhood maltreatment is an important risk factor for adult depression and has been associated with changes in the hypothalamic pituitary adrenal (HPA) axis, including cortisol secretion and methylation of the FKBP5 gene. Furthermore, associations between depression and HPA changes have been reported. This study investigated the associations of whole-blood FKBP5 mRNA levels, serum cortisol levels, childhood maltreatment, and depressive symptoms with the whole-blood methylation status (assessed via target bisulfite sequencing) of 105 CpGs at the FKBP5 locus using data from the general population-based Study of Health in Pomerania (SHIP) (N = 203). Both direct and interaction effects with the rs1360780 single-nucleotide polymorphism were investigated. Nominally significant associations of main effects on methylation of a single CpG site were observed at intron 3, intron 7, and the 3'-end of the gene. Additionally, methylation at two clusters at the 3'-end and intron 7 were nominally associated with childhood maltreatment × rs1360780 and depressive symptoms × rs1360780, respectively. The results add to the understanding of molecular mechanisms underlying the emergence of depression and could aid the development of personalised depression therapy and drug development.

Polygenic risk for schizophrenia converges on alternative polyadenylation as molecular mechanism underlying synaptic impairment.

Schizophrenia (SCZ) is a genetically heterogenous psychiatric disorder of highly polygenic nature. Correlative evidence from genetic studies indicate that the aggregated effects of distinct genetic risk factor combinations found in each patient converge onto common molecular mechanisms. To prove this on a functional level, we employed a reductionistic cellular model system for polygenic risk by differentiating induced pluripotent stem cells (iPSCs) from 104 individuals with high polygenic risk load and controls into cortical glutamatergic neurons (iNs). Multi-omics profiling identified widespread differences in alternative polyadenylation (APA) in the 3' untranslated region of many synaptic transcripts between iNs from SCZ patients and healthy donors. On the cellular level, 3'APA was associated with a reduction in synaptic density of iNs. Importantly, differential APA was largely conserved between postmortem human prefrontal cortex from SCZ patients and healthy donors, and strongly enriched for transcripts related to synapse biology. 3'APA was highly correlated with SCZ polygenic risk and affected genes were significantly enriched for SCZ associated common genetic variation. Integrative functional genomic analysis identified the RNA binding protein and SCZ GWAS risk gene PTBP2 as a critical trans-acting factor mediating 3'APA of synaptic genes in SCZ subjects. Functional characterization of PTBP2 in iNs confirmed its key role in 3'APA of synaptic transcripts and regulation of synapse density. Jointly, our findings show that the aggregated effects of polygenic risk converge on 3'APA as one common molecular mechanism that underlies synaptic impairments in SCZ.

Extensive evaluation of DNA methylation of functional elements in the murine Fkbp5 locus using high-accuracy DNA methylation measurement via targeted bisulfite sequencing.

FKBP5 is an important stress-regulatory gene implicated in stress-related psychiatric diseases. Single nucleotide polymorphisms of the FKBP5 gene were shown to interact with early life stress to alter the glucocorticoid-related stress response and moderate disease risk. Demethylation of cytosine-phosphate-guanine-dinucleotides (CpGs) in regulatory glucocorticoid-responsive elements was suggested to be the mediating epigenetic mechanism for long-term stress effects, but studies on Fkbp5 DNA methylation (DNAm) in rodents are so far limited. We evaluated the applicability of high-accuracy DNA methylation measurement via targeted bisulfite sequencing (HAM-TBS), a next-generation sequencing-based technology, to allow a more in-depth characterisation of the DNA methylation of the murine Fkbp5 locus in three different tissues (blood, frontal cortex and hippocampus). In this study, we not only increased the number of evaluated sites in previously described regulatory regions (in introns 1 and 5), but also extended the evaluation to novel, possibly relevant regulatory regions of the gene (in intron 8, the transcriptional start site, the proximal enhancer and CTCF-binding sites within the 5'UTR). We here describe the assessment of HAM-TBS assays for a panel of 157 CpGs with possible functional relevance in the murine Fkbp5 gene. DNAm profiles were tissue-specific, with lesser differences between the two brain regions than between the brain and blood. Moreover, we identified DNAm changes in the Fkbp5 locus after early life stress exposure in the frontal cortex and blood. Our findings indicate that HAM-TBS is a valuable tool for broader exploration of the DNAm of the murine Fkbp5 locus and its involvement in the stress response.

Transdiagnostic evaluation of epigenetic age acceleration and burden of psychiatric disorders.

Different psychiatric disorders as well as exposure to adverse life events have individually been associated with multiple age-related diseases and mortality. Age acceleration in different epigenetic clocks can serve as biomarker for such risk and could help to disentangle the interplay of psychiatric comorbidity and early adversity on age-related diseases and mortality. We evaluated five epigenetic clocks (Horvath, Hannum, PhenoAge, GrimAge and DunedinPoAm) in a transdiagnostic psychiatric sample using epigenome-wide DNA methylation data from peripheral blood of 429 subjects from two studies at the Max Planck Institute of Psychiatry. Burden of psychiatric disease, represented by a weighted score, was significantly associated with biological age acceleration as measured by GrimAge and DunedinPoAm (R2-adj. 0.22 and 0.33 for GrimAge and DunedinPoAm, respectively), but not the other investigated clocks. The relation of burden of psychiatric disease appeared independent of differences in socioeconomic status and medication. Our findings indicate that increased burden of psychiatric disease may associate with accelerated biological aging. This highlights the importance of medical management of patients with multiple psychiatric comorbidities and the potential usefulness of specific epigenetic clocks for early detection of risk and targeted intervention to reduce mortality in psychiatric patients.

DiffBrainNet: Differential analyses add new insights into the response to glucocorticoids at the level of genes, networks and brain regions.

Genome-wide gene expression analyses are invaluable tools for studying biological and disease processes, allowing a hypothesis-free comparison of expression profiles. Traditionally, transcriptomic analysis has focused on gene-level effects found by differential expression. In recent years, network analysis has emerged as an important additional level of investigation, providing information on molecular connectivity, especially for diseases associated with a large number of linked effects of smaller magnitude, like neuropsychiatric disorders. Here, we describe how combined differential expression and prior-knowledge-based differential network analysis can be used to explore complex datasets. As an example, we analyze the transcriptional responses following administration of the glucocorticoid/stress receptor agonist dexamethasone in 8 mouse brain regions important for stress processing. By applying a combination of differential network- and expression-analyses, we find that these explain distinct but complementary biological mechanisms of the glucocorticoid responses. Additionally, network analysis identifies new differentially connected partners of risk genes and can be used to generate hypotheses on molecular pathways affected. With DiffBrainNet (http://diffbrainnet.psych.mpg.de), we provide an analysis framework and a publicly available resource for the study of the transcriptional landscape of the mouse brain which can identify molecular pathways important for basic functioning and response to glucocorticoids in a brain-region specific manner.

FKBP5 intron 7 methylation is associated with higher anxiety proneness and smaller right thalamus volume in adolescents.

Dysregulation of stress response systems may mediate the detrimental effects of childhood trauma (CT) on mental health. FKBP5 regulates glucocorticoid receptor sensitivity and exerts pleiotropic effects on intracellular signaling, neurobiology and behavior. We investigated whether CT, alone and in combination with rs1360780 genotype, is associated with altered FKBP5 methylation and whether CT-associated methylation profiles are associated with anxiety proneness (AP) and structural brain volumes. Ninety-four adolescents completed the Childhood Trauma Questionnaire, and a composite AP score was generated from the Childhood Anxiety Sensitivity Index and the State-Trait Anxiety Inventory-Trait measure. Mean methylation values for 12 regulatory regions and 25 individual CpG sites were determined using high-accuracy measurement via targeted bisulfite sequencing. FKBP5 rs1360780 genotype and structural MRI data were available for a subset of participants (n = 71 and n = 75, respectively). Regression models revealed an inverse association between methylation of three intron 7 CpG sites (35558438, 35558566 and 35558710) and right thalamus volume. CpG35558438 methylation was positively associated with AP scores. Our data indicate that an intron 7 methylation profile, consistent with lower FKBP5 expression and elevated high sensitivity glucocorticoid receptor levels, is associated with higher AP and smaller right thalamus volume. Research into the mechanisms underlying the intron 7 methylation-thalamus volume relationship, and whether it confers increased risk for long-term psychopathology by altering the regulatory threshold of stress responding, is required.

Genome-Wide Copy Number Variant and High-Throughput Transcriptomics Analyses of Placental Tissues Underscore Persisting Child Susceptibility in At-Risk Pregnancies Cleared in Standard Genetic Testing.

Several studies have shown that children from pregnancies with estimated first-trimester risk based on fetal nuchal translucency thickness and abnormal maternal serum pregnancy protein and hormone levels maintain a higher likelihood of adverse outcomes, even if initial testing for known genetic conditions is negative. We used the Finnish InTraUterine cohort (ITU), which is a comprehensively characterized perinatal cohort consisting of 943 mothers and their babies followed throughout pregnancy and 18 months postnatally, including mothers shortlisted for prenatal genetic testing but cleared for major aneuploidies (cases: n = 544, 57.7%) and control pregnancies (n = 399, 42.3%). Using genome-wide genotyping and RNA sequencing of first-trimester and term placental tissue, combined with medical information from registry data and maternal self-report data, we investigated potential negative medical outcomes and genetic susceptibility to disease and their correlates in placenta gene expression. Case mothers did not present with higher levels of depression, perceived stress, or anxiety during pregnancy. Case children were significantly diagnosed more often with congenital malformations of the circulatory system (4.12 (95% CI [1.22−13.93]) higher hazard) and presented with significantly more copy number duplications as compared to controls (burden analysis, based on all copy number variants (CNVs) with at most 10% frequency, 823 called duplications in 297 cases versus 626 called duplications in 277 controls, p = 0.01). Fifteen genes showed differential gene expression (FDR < 0.1) in association with congenital malformations in first-trimester but not term placenta. These were significantly enriched for genes associated with placental dysfunction. In spite of normal routine follow-up prenatal testing results in early pregnancy, case children presented with an increased likelihood of negative outcomes, which should prompt vigilance in follow-up during pregnancy and after birth.

Cohort profile: InTraUterine sampling in early pregnancy (ITU), a prospective pregnancy cohort study in Finland: study design and baseline characteristics.

PURPOSE: The InTraUterine sampling in early pregnancy (ITU) is a prospective pregnancy cohort study. The overarching aim of ITU is to unravel genomic, epigenomic, transcriptomic, endocrine, inflammatory and metabolic maternal-placental-fetal mechanisms involved in the programming of health and disease after exposure to prenatal environmental adversity, such as maternal malnutrition, cardiometabolic disorders, infections, medical interventions, mental disorders and psychosocial stress. This paper describes the study protocol, design and baseline characteristics of the cohort. PARTICIPANTS: We included 944 pregnant Finnish women, their partners and children born alive between April 2012 and December 2017. The women were recruited through the national, voluntary trisomy 21 screening between 9+0 and 21+6 gestational weeks. Of the participating women, 543 were screen positive and underwent fetal chromosomal testing. Test result of these women suggested no fetal chromosomal abnormality. Further, we recruited 401 women who were screen negative and who did not undergo fetal chromosomal testing. FINDINGS TO DATE: We have collected chorionic villi and amniotic fluid from the screen-positive women; blood, urine, buccal swabs and diurnal salivary samples from all women; blood and buccal swabs from all partners; and placenta, cord blood and buccal swabs from all newborns for analyses of the genome, epigenome, transcriptome, and endocrine, inflammatory and metabolic markers. These data are coupled with comprehensive phenotypes, including questions on demographic characteristics, health and well-being of the women and their partners during pregnancy and of the women and their children at the child's age of 1.7 and 3 years. Data also come from patient records and nationwide registers covering health, lifestyle and medication data. FUTURE PLANS: Multiple layers of ITU data allow integrative data analyses, which translate to biomarker identification and allow risk stratification and understanding of the biological mechanisms involved in prenatal programming of health and disease.

Cell-Type-Specific Impact of Glucocorticoid Receptor Activation on the Developing Brain: A Cerebral Organoid Study.

OBJECTIVE: A fine-tuned balance of glucocorticoid receptor (GR) activation is essential for organ formation, with disturbances influencing many health outcomes. In utero, glucocorticoids have been linked to brain-related negative outcomes, with unclear underlying mechanisms, especially regarding cell-type-specific effects. An in vitro model of fetal human brain development, induced human pluripotent stem cell (hiPSC)-derived cerebral organoids, was used to test whether cerebral organoids are suitable for studying the impact of prenatal glucocorticoid exposure on the developing brain. METHODS: The GR was activated with the synthetic glucocorticoid dexamethasone, and the effects were mapped using single-cell transcriptomics across development. RESULTS: The GR was expressed in all cell types, with increasing expression levels through development. Not only did its activation elicit translocation to the nucleus and the expected effects on known GR-regulated pathways, but also neurons and progenitor cells showed targeted regulation of differentiation- and maturation-related transcripts. Uniquely in neurons, differentially expressed transcripts were significantly enriched for genes associated with behavior-related phenotypes and disorders. This human neuronal glucocorticoid response profile was validated across organoids from three independent hiPSC lines reprogrammed from different source tissues from both male and female donors. CONCLUSIONS: These findings suggest that excessive glucocorticoid exposure could interfere with neuronal maturation in utero, leading to increased disease susceptibility through neurodevelopmental processes at the interface of genetic susceptibility and environmental exposure. Cerebral organoids are a valuable translational resource for exploring the effects of glucocorticoids on early human brain development.

The pediatric buccal epigenetic clock identifies significant ageing acceleration in children with internalizing disorder and maltreatment exposure.

BACKGROUND: Studies reporting accelerated ageing in children with affective disorders or maltreatment exposure have relied on algorithms for estimating epigenetic age derived from adult samples. These algorithms have limited validity for epigenetic age estimation during early development. We here use a pediatric buccal epigenetic (PedBE) clock to predict DNA methylation-based ageing deviation in children with and without internalizing disorder and assess the moderating effect of maltreatment exposure. We further conduct a gene set enrichment analysis to assess the contribution of glucocorticoid signaling to PedBE clock-based results. METHOD: DNA was isolated from saliva of 158 children [73 girls, 85 boys; mean age (SD) = 4.25 (0.8) years] including children with internalizing disorder and maltreatment exposure. Epigenetic age was estimated based on DNA methylation across 94 CpGs of the PedBE clock. Residuals of epigenetic age regressed against chronological age were contrasted between children with and without internalizing disorder. Maltreatment was coded in 3 severity levels and entered in a moderation model. Genome-wide dexamethasone-responsive CpGs were derived from an independent sample and enrichment of these CpGs within the PedBE clock was identified. RESULTS: Children with internalizing disorder exhibited significant acceleration of epigenetic ageing as compared to children without internalizing disorder (F1,147 = 6.67, p = .011). This association was significantly moderated by maltreatment severity (b = 0.49, 95% CI [0.073, 0.909], t = 2.322, p = .022). Children with internalizing disorder who had experienced maltreatment exhibited ageing acceleration relative to children with no internalizing disorder (1-2 categories: b = 0.50, 95% CI [0.170, 0.821], t = 3.008, p = .003; 3 or more categories: b = 0.99, 95% CI [0.380, 1.593], t = 3.215, p = .002). Children with internalizing disorder who were not exposed to maltreatment did not show epigenetic ageing acceleration. There was significant enrichment of dexamethasone-responsive CpGs within the PedBE clock (OR = 4.36, p = 1.65*10-6). Among the 94 CpGs of the PedBE clock, 18 (19%) were responsive to dexamethasone. CONCLUSION: Using the novel PedBE clock, we show that internalizing disorder is associated with accelerated epigenetic ageing in early childhood. This association is moderated by maltreatment severity and may, in part, be driven by glucocorticoids. Identifying developmental drivers of accelerated epigenetic ageing after maltreatment will be critical to devise early targeted interventions.

Childhood adversity correlates with stable changes in DNA methylation trajectories in children and converges with epigenetic signatures of prenatal stress.

Childhood maltreatment (CM) is an established major risk factor for a number of negative health outcomes later in life. While epigenetic mechanisms, such as DNA methylation (DNAm), have been proposed as a means of embedding this environmental risk factor, little is known about its timing and trajectory, especially in very young children. It is also not clear whether additional environmental adversities, often experienced by these children, converge on similar DNAm changes. Here, we calculated a cumulative adversity score, which additionally to CM includes socioeconomic status (SES), other life events, parental psychopathology and epigenetic biomarkers of prenatal smoking and alcohol consumption. We investigated the effects of CM alone as well as the adversity score on longitudinal DNAm trajectories in the Berlin Longitudinal Child Study. This is a cohort of 173 children aged 3-5 years at baseline of whom 86 were exposed to CM. These children were followed-up for 2 years with extensive psychometric and biological assessments as well as saliva collection at 5 time points providing genome-wide DNAm levels. Overall, only a few DNAm patterns were stable over this timeframe, but less than 10 DNAm regions showed significant changes. At baseline, neither CM nor the adversity score associated with DNAm changes. However, in 6 differentially methylated regions (DMRs), CM and the adversity score significantly moderated DNAm trajectories over time. A number of these DMRs have previously been associated with adverse prenatal exposures. In our study, children exposed to CM also presented with epigenetic signatures indicative of increased prenatal exposure to tobacco and alcohol, as compared to non-CM exposed children. These epigenetic signatures of prenatal exposure strongly correlate with DNAm regions associated with CM and the adversity score. Finally, weighted correlation network analysis revealed a module of CpGs exclusively associated with CM. While our study identifies DNAm loci specifically associated with CM, especially within long non-coding RNAs, the majority of associations were found with the adversity score with convergent association with indicators of adverse prenatal exposures. This study highlights the importance of mapping not only of the epigenome but also the exposome and extending the observational timeframe to well before birth.

The Atherosclerosis Risk Variant rs2107595 Mediates Allele-Specific Transcriptional Regulation of HDAC9 via E2F3 and Rb1.

Background and Purpose- Genome-wide association studies have identified the HDAC9 (histone deacetylase 9) gene region as a major risk locus for atherosclerotic stroke and coronary artery disease in humans. Previous results suggest a role of altered HDAC9 expression levels as the underlying disease mechanism. rs2107595, the lead single nucleotide polymorphism for stroke and coronary artery disease resides in noncoding DNA and colocalizes with histone modification marks suggestive of enhancer elements. Methods- To determine the mechanisms by which genetic variation at rs2107595 regulates HDAC9 expression and thus vascular risk we employed targeted resequencing, proteome-wide search for allele-specific nuclear binding partners, chromatin immunoprecipitation, genome-editing, reporter assays, circularized chromosome conformation capture, and gain- and loss-of-function experiments in cultured human cell lines and primary immune cells. Results- Targeted resequencing of the HDAC9 locus in patients with atherosclerotic stroke and controls supported candidacy of rs2107595 as the causative single nucleotide polymorphism. A proteomic search for nuclear binding partners revealed preferential binding of the E2F3/TFDP1/Rb1 complex (E2F transcription factor 3/transcription factor Dp-1/Retinoblastoma 1) to the rs2107595 common allele, consistent with the disruption of an E2F3 consensus site by the risk allele. Gain- and loss-of-function studies showed a regulatory effect of E2F/Rb proteins on HDAC9 expression. Compared with the common allele, the rs2107595 risk allele exhibited higher transcriptional capacity in luciferase assays and was associated with higher HDAC9 mRNA levels in primary macrophages and genome-edited Jurkat cells. Circularized chromosome conformation capture revealed a genomic interaction of the rs2107595 region with the HDAC9 promoter, which was stronger for the common allele as was the in vivo interaction with E2F3 and Rb1 determined by chromatin immunoprecipitation. Gain-of-function experiments in isogenic Jurkat cells demonstrated a key role of E2F3 in mediating rs2107595-dependent transcriptional regulation of HDAC9. Conclusions- Collectively, our findings imply allele-specific transcriptional regulation of HDAC9 via E2F3 and Rb1 as a major mechanism mediating vascular risk at rs2107595.

Polygenic Risk: Predicting Depression Outcomes in Clinical and Epidemiological Cohorts of Youths.

OBJECTIVE: Identifying risk factors for major depression and depressive symptoms in youths could have important implications for prevention efforts. This study examined the association of polygenic risk scores (PRSs) for a broad depression phenotype derived from a large-scale genome-wide association study (GWAS) in adults, and its interaction with childhood abuse, with clinically relevant depression outcomes in clinical and epidemiological youth cohorts. METHODS: The clinical cohort comprised 279 youths with major depression (mean age=14.76 years [SD=2.00], 68% female) and 187 healthy control subjects (mean age=14.67 years [SD=2.45], 63% female). The first epidemiological cohort included 1,450 youths (mean age=13.99 years [SD=0.92], 63% female). Of those, 694 who were not clinically depressed at baseline underwent follow-ups at 6, 12, and 24 months. The replication epidemiological cohort comprised children assessed at ages 8 (N=184; 49.2% female) and 11 (N=317; 46.7% female) years. All cohorts were genome-wide genotyped and completed measures for major depression, depressive symptoms, and/or childhood abuse. Summary statistics from the largest GWAS to date on depression were used to calculate the depression PRS. RESULTS: In the clinical cohort, the depression PRS predicted case-control status (odds ratio=1.560, 95% CI=1.230-1.980), depression severity (ß=0.177, SE=0.069), and age at onset (ß=-0.375, SE=0.160). In the first epidemiological cohort, the depression PRS predicted baseline depressive symptoms (ß=0.557, SE=0.200) and prospectively predicted onset of moderate to severe depressive symptoms (hazard ratio=1.202, 95% CI=1.045-1.383). The associations with depressive symptoms were replicated in the second epidemiological cohort. Evidence was found for an additive, but not an interactive, effect of the depression PRS and childhood abuse on depression outcomes. CONCLUSIONS: Depression PRSs derived from adults generalize to depression outcomes in youths and may serve as an early indicator of clinically significant levels of depression.

Evidence for oestrogen sensitivity in perinatal depression: pharmacological sex hormone manipulation study.

BACKGROUND: Enhanced sensitivity to oestrogen signalling may drive increased risk for depressive symptoms when exposed to peripartum sex-steroid hormone fluctuations. AIM: Testing if 116 pre-identified sex steroid-responsive transcripts that predicted perinatal depression (PND) translates to a pharmacological model of hormone-induced mood changes. METHOD: We generated longitudinal, genome-wide gene-expression and DNA-methylation data from 60 women exposed to a gonadotrophin-releasing hormone agonist (GnRHa) or placebo. We used linear mixed-effect models to assess differences between baseline and follow-up for gene expression and DNA methylation in the biphasic ovarian response to GnRHa. RESULTS: Of the 116 PND-predictive transcripts, a significant (19%) overlap was observed with those differentially expressed post-GnRHa at both early and later follow-up, indicating sustained effects. Similarly, 49% of tested genes were differentially methylated post-GnRHa at the late follow-up. Within the GnRHa group, a large proportion of PND genes were significantly associated (gene expression; DNA methylation) with changes in depressive symptoms (28%; 66%), oestradiol levels (49%; 66%) and neocortex serotonin transporter binding (8%; 45%) between baseline and follow-up. CONCLUSIONS: Our data bridge clinical PND biomarkers with a pharmacological model of sex hormone-induced mood changes and directly relate oestrogen-induced biological changes with depressive symptoms and associated serotonin-signalling changes. Our data highlight that individual variations in molecular sensitivity to oestrogen associate with susceptibility to hormone-induced mood changes and hold promise for candidate biomarkers. DECLARATION OF INTEREST: V.G.F. received honorarium for being a speaker for H. Lundbeck A/S. E.B.B. receives research funding from Böhringer Ingelheim to investigate FKBP5 as a potential drug target for depression.

The Role of m6A/m-RNA Methylation in Stress Response Regulation.

N6-methyladenosine (m6A) and N6,2'-O-dimethyladenosine (m6Am) are abundant mRNA modifications that regulate transcript processing and translation. The role of both, here termed m6A/m, in the stress response in the adult brain in vivo is currently unknown. Here, we provide a detailed analysis of the stress epitranscriptome using m6A/m-seq, global and gene-specific m6A/m measurements. We show that stress exposure and glucocorticoids region and time specifically alter m6A/m and its regulatory network. We demonstrate that deletion of the methyltransferase Mettl3 or the demethylase Fto in adult neurons alters the m6A/m epitranscriptome, increases fear memory, and changes the transcriptome response to fear and synaptic plasticity. Moreover, we report that regulation of m6A/m is impaired in major depressive disorder patients following glucocorticoid stimulation. Our findings indicate that brain m6A/m represents a novel layer of complexity in gene expression regulation after stress and that dysregulation of the m6A/m response may contribute to the pathophysiology of stress-related psychiatric disorders.

Sex-related differential response to dexamethasone in endocrine and immune measures in depressed in-patients and healthy controls.

Although sex differences in major depression have been reported repeatedly, the underlying mechanisms are still disputed. The rapidly changing gonadal steroid concentrations of the postpartum period or during menopause have been shown to be associated with depressive symptoms and to modulate the hypothalamic-pituitary-adrenal (HPA)-axis, which is implicated in depression. The sample comprised of 128 depressed in-patients (36.7% women) and 166 healthy controls (30.0% women). Blood was collected at baseline (at 6pm) and then 3 h as well as 21 h after ingestion of 1.5 mg dexamethasone for measurement of cortisol, ACTH and blood count. To further assess the function of the HPA-axis the dexamethasone/corticotrophin releasing hormone (Dex-CRH) test was performed in a subsample of 115 patients and 116 controls the following day. A significant interaction effect between sex, disease and ACTH concentrations over time after dexamethasone stimulation was observed, with men showing increased ACTH concentrations at baseline and after 21 h, while there was no difference after 3 h (p = .007). After separating for disease status this significant interaction effect was only observed in controls (p = .005). The cortisol response in the dex-CRH test was enhanced in female compared to male controls (p = .002). Leucocytes showed a stronger increase upon dexamethasone administration only in female compared to male controls (p = .023). These findings suggest a higher glucocorticoid receptor sensitivity following in-vivo glucocorticoid stimulation in healthy women that was absent in depressed patients. The sex-related differences in HPA-axis regulation and immune system function may contribute to the vulnerability of female sex to the development of depression.

Common genes associated with antidepressant response in mouse and man identify key role of glucocorticoid receptor sensitivity.

Response to antidepressant treatment in major depressive disorder (MDD) cannot be predicted currently, leading to uncertainty in medication selection, increasing costs, and prolonged suffering for many patients. Despite tremendous efforts in identifying response-associated genes in large genome-wide association studies, the results have been fairly modest, underlining the need to establish conceptually novel strategies. For the identification of transcriptome signatures that can distinguish between treatment responders and nonresponders, we herein submit a novel animal experimental approach focusing on extreme phenotypes. We utilized the large variance in response to antidepressant treatment occurring in DBA/2J mice, enabling sample stratification into subpopulations of good and poor treatment responders to delineate response-associated signature transcript profiles in peripheral blood samples. As a proof of concept, we translated our murine data to the transcriptome data of a clinically relevant human cohort. A cluster of 259 differentially regulated genes was identified when peripheral transcriptome profiles of good and poor treatment responders were compared in the murine model. Differences in expression profiles from baseline to week 12 of the human orthologues selected on the basis of the murine transcript signature allowed prediction of response status with an accuracy of 76% in the patient population. Finally, we show that glucocorticoid receptor (GR)-regulated genes are significantly enriched in this cluster of antidepressant-response genes. Our findings point to the involvement of GR sensitivity as a potential key mechanism shaping response to antidepressant treatment and support the hypothesis that antidepressants could stimulate resilience-promoting molecular mechanisms. Our data highlight the suitability of an appropriate animal experimental approach for the discovery of treatment response-associated pathways across species.

Sequencing on the SOLiD 5500xl System - in-depth characterization of the GC bias.

Different types of sequencing biases have been described and subsequently improved for a variety of sequencing systems, mostly focusing on the widely used Illumina systems. Similar studies are missing for the SOLiD 5500xl system, a sequencer which produced many data sets available to researchers today. Describing and understanding the bias is important to accurately interpret and integrate these published data in various ongoing research projects. We report a particularly strong GC bias for this sequencing system when analyzing a defined gDNA mix of 5 microbes with a wide range of different GC contents (20-72%) when comparing to the expected distribution and Illumina MiSeq data from the same DNA pool. Since we observed this bias already under PCR-free conditions, changing the PCR conditions during library preparation - a common strategy to handle bias in the Illumina system - was not relevant. Source of the bias appeared to be an uneven heat distribution during the SOLiD emulsion PCR (ePCR) - for enrichment of libraries prior loading - since ePCR in either small pouches or in 96-well plates improved the GC bias. Sequencing of chromatin immunoprecipitated DNA (ChIP-seq) is a common approach in epigenetics. ChIP-seq of the mixed source histone mark H3K9ac (acetyl Histone H3 lysine 9), typically found on promoter regions and on gene bodies, including CpG islands, performed on a SOLiD 5500xl machine, resulted in major loss of reads at GC rich loci (GC content ≥ 62%), not explained by low sequencing depth. This was improved with adaptations of the ePCR.